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^ 1 Abstract 
(N ; 

We consider a magnetic impurity which interacts by hybridization with a sys- 

■ tern of strongly correlated conduction electrons. The latter are described by 
CN ■ a Hubbard Hamiltonian. By means of a canconical transformation the charge 

o 

■ degrees of freedom of the magnetic impurity are eliminated. The resulting 
effective Hamiltonian H e g is investigated and various limiting cases are con- 

a 

g 

o 

o 



sidered. If the Hubbard interaction U between the conduction electrons is 
neglected H e R reduces to a form obtained by the Schrieffer- Wolff transfor- 
mation, which is essentially the Kondo Hamiltonian. If U is large and the 



^ . correlations are strong H e fi is changed. One modification concerns the coef- 

ficient of the dominant exchange coupling of the magnetic impurity with the 
nearest lattice site. When the system is hole doped, there is also an antiferro- 
magnetic coupling to the nearest neighbors of that site involving additionally 
a hole. Furthermore, it is found that the magnetic impurity attracts a hole. 
In the case of electron doping, double occupancies are repelled by the impu- 
rity. In contrast to the hole-doped case, we find no magnetic coupling which 

additionally involves a doubly occupied site. 
71.27.+a, 75.30.Hx, 71.70.Gm 
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I. INTRODUCTION 



Recently, heavy-fermion behavior has been observed in the electron-doped cuprate 
Nd^Ce^CuC^ (0.1 < x < 0.2)1 Below 0.3 K a linear specific heat C v = 7T is observed 
with a large Sommerfeld coefficient 7 ~ 4 J/ (mole Nd • K 2 ). In the same temperature regime, 
the spin susceptibility is found to be independent of the temperature and the Sommerfeld- 
Wilson ratio is of order unity. These are characteristic features of heavy-fermion excitations.ll 
However, the characteristic low energy scale of the order of 1 K which is associated with 
this behavior is not based on a Kondo lattice effect, as it is the case in other heavy-fermion 
systems. Rather, it is based on a Zeeman effect, a consequence of the strong electron corre- 
lations in the conducting Cu0 2 planes.! Indeed, undoped Nd 2 Cu0 4 is an antiferromagnetic 
charge-transfer insulator instead of a metal,B! despite of one hole per unit cell. The Nd ions 
are therefore coupling to a system of strongly correlated electrons^ rather than to weakly 
or uncorrelated ones. 

As a first step towards an understanding of the consequences of the strong correlations, 
the effect of a single impurity should be investigated. Several authors studied the influence of 
non-magnetic impurities in systems with strongly correlated electron in order to explain the 
magnetic properties of the undoped or weakly doped host materials of the high-temperature 
superconductors (e.g., La 2 Cu0 4 ). The considered impurities stem from substituting Sr 
for La,iF0 or from substitutions in the Cu0 2 planes, i.e., Zn for CuEOHlI and S for 0.0 
Naeaosa et aM treated magnetic impurities in an undoped two-dimensional Heisenberg 
antiferromagnet by adding external spins with a different exchange coupling constant. Nagao 
et a/.lll discussed the influence of both, non-magnetic and magnetic impurity scattering on 
the spin density wave state and the superconducting state within a Hubbard model for 
weakly correlated conduction electrons. 

In this paper, we study the coupling of a single magnetic impurity to strongly correlated 
electrons moving on a lattice. The latter are described by a Hubbard Hamiltonian and the 
total Hamiltonian must, therefore, go beyond that of the single-site Anderson impurity 
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model.O By means of a canonical transformation the charge degrees of freedom of the 
impurity site are eliminated like by the Schrieffer- Wolff transformational in the case of 
uncorrelated electrons. Due to the strong correlations of the conduction electrons, new 
terms appear in the resulting effective Hamiltonian. They are analyzed and interpreted. 
Special attention is paid to the case of a nearly half-filled conduction band with strong 
correlations, a situation prevailing, e.g., in Nd2-zCe a ;CuC)4. We believe that the new terms 
resulting from the correlations of the conduction electrons are of relevance for a number of 
different systems. If the correlations of the conduction electrons are neglected the effective 
Hamiltonian reduces again to that of Schrieffer and Wolff.El 

In the next section the starting Hamiltonian is defined and the canonical transformation 
is outlined. The resulting effective Hamiltonian for uncorrelated electrons is discussed in 
Sec. |TJ. Particularly interesting is the limit of strong correlations. In Sec. |TV| , we discuss 
the system at half filling and both, the hole and the electron doped case. Finally, Sec. [V| 
contains a summary of the results and the conclusions. 

II. THE HAMILTONIAN AND ITS TRANSFORMATION 

We describe the strongly correlated conduction electrons in the substrate by a Hubbard 
Hamiltonian on a hypercubic lattice with unit vectors x 

H c = -tJ2 + U £ n n n rl = H t + H v . (2.1) 

jx,cr j 

The operators c] CT (c Jcr ) create (destroy) an electron with spin a on site j and rij a = c^ a Cj a . 
The magnetic impurity is assumed to contain one orbital (e.g., 4/), which is either empty 
or singly occupied. Double occupancies are excluded because of the strong repulsion of 
electrons in that orbital. The orbital energy is therefore given by 

Hf = e f Y,ith (2.2) 

with fl = fl{l — rif), rif = J2a f\fai an d e/ < 0. The operators f\ and f a obey the usual 
anticommutation relations. For later convenience, we define r = —ej/U(> 0). Finally, we 
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assume that the interaction between the / orbital and the conduction electrons is local and 
described by a hybridization contribution like in the Anderson mo del§ 

^-/^E^i + Zk)- (2-3) 

Instead of coupling the / orbital with the site of the Hubbard lattice we could have replaced 
one site of the lattice by the / orbital. The present choice is motivated by the situation 
prevailing in Nd2CuC>4. 

The total Hamiltonian reads 

H — H c + Hf + H c -f 

= - t Y, C ^ C j+xa + UY, n ^ n 3l + e fY,flL (2.4) 

jx,cr j cr 

a 

From this Hamiltonian we want to derive an effective one, H e s, which acts on the space 
with a singly occupied / orbital, only. Configurations with an empty / orbital are elim- 
inated by first performing a canonical transformation, thereby eliminating the mixing of 
configurations with different / orbital occupancies and projecting afterwards onto the space 
with Uj = 1. The canonical transformation is written as 

# can = e s He~ s , (2.5) 

where S is determined by the requirement that H c _f disappears to lowest order in V. This 
leads to 

[H c + H f , S\. = (L c + L f )S = H c _ f . (2.6) 

Here we introduced the Liouville operators L c and Lf which act on an operator A according 
to L C A = [H c , A]_ and LfA = [Hf, A]_, respectively. Equation (p.6|) has the formal solution 

and up to terms of second order in V the Hamiltonian H can reads 



Next, we project by means of a projector Pf onto the space with n/ = 1. The effective 
Hamiltonian is then given by 

-ffcff — PfHc&nPf 

= e/ + H c + — E(/a'/ CT c 0ff , -cl + he). 

<t,<x J L 

To proceed further we use the expansion 
1 1 



51 ( : — i~ Lt ) - — j— c Li 



(2.10) 



which is one in powers of t/ef, and terminate it after v — 2. The Liouvilleans L t and Ljj in 
Eq. ( |2.10| ) are defined by L t A = [H t ,A]- and LyA = [Hu,A}_, respectively. These terms 
are easily evaluated if one expresses c| CT in terms of eigenoperators of Lu, i.e., 

Ac = Ad 1 - ms) + c \a n is = Ac + Ac (2-11) 

where a = —a. The last two operators have the eigenvalues and 1, respectively. The 
resulting terms of Eq. (|2.10|) up to second order are listed in Appendix |A], where next- 
nearest neighbor contributions which arise in second order [y = 2) have been neglected. 
Inserting these terms into Eq. ( [2.9|) one obtains if e ff in the form 

H cS = e f + H c + + + H W . (2.12) 

The Hamiltonians are of order t v . The general expressions, which hold for arbitrary 
values of U/ef, are lengthy and are, therefore, moved to Appendix [A[ In the following we 
discuss and list only special limits of them. 

III. LIMIT OF UNCORRELATED CONDUCTION ELECTRONS 



For U = (r = oo), we obtain from Eq. ( |2.12| ) the result of the conventional Schrieffer- 
Wolff transformation@, i.e., 
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y\/ 2 

H eS (U = 0) = e f +<—(2-n ) 

2e f 



~l~ ( C jcr C j+XCT "I" C j+XCT C jo-) (3-1) 

27V 2 t\/ 2 

S/s H — S/ XX Si ° + s o*)- 

e / e / i(0) 

Here, the summation z(0) is over the Z nearest neighbor sites of site 0, Sf is the / electron 
spin, s is that of site and 

Z a/3 

cr denote the Pauli matrices. 

The last two terms in Eq. ( |3.1| ) describe the Kondo coupling between the / spin and 
the spin of the conduction electrons.^ Applying L t in Eq. ( |2.10| ) results in a delocalization 
of the interaction to first order in t. To second order, we obtain a renormalization factor 
7' = 1 + Zt 2 /e"j for S/So- Further delocalization of the interaction is neglected in our 
approximation, since next-nearest neighbor contributions were not taken into account in the 
truncated form of Eq. (|2.10|) , see Appendix [A]. Due to the impurity the electron hopping is 
reduced between site and its nearest neighbors 

t (_ V*_ 

Finally, there is a one-particle potential, which reduces the probability of having an electron 
at site 0, as this would diminish the gain in kinetic energy of the / electron by a virtual hop. 



^° ) = -2( 1 -^o|. (3.3) 



IV. LIMIT OF STRONG CORRELATIONS 

Of particular interest is the limit of strong correlations (U ^> t), which we treat in the 
following. In contrast to the ordinary Hubbard model (without impurity), the two cases 
of doping the half-filled system with holes or with electrons are no longer symmetric since 
we take into account an empty and a singly occupied / orbital only. In the following, we 
perform a degenerate perturbation theory for both cases and comment on the half-filled 
case. 
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A. Hole doping 



For less than half filling (hole doping), doubly occupied sites are excluded in the limit 
of strong correlations (large U). Let P c denote the projector onto the configuration space 
without doubly occupied sites. Furthermore, set Q c = 1 — P c . The effective Hamiltonian, 
when reduced to the space of no doubly occupied lattice sites become s@ 

# (h) = P c H cS P c + P c H cS Q c Q c H cS P c (4.1) 

where is the ground-state energy. Only those parts of H e Q contribute to the last term 
to order U^ 1 which generate a single doubly occupied site. Thus, 

# (h) = P c H cS P c - ^Pc(H t + H^)Q c (H t + H^)P C (4.2) 

(see Appendix [AJ). When terms of order V 4 , as well as terms involving three lattice sites are 
neglected the second term in Eq. ( |4.2|) becomes 



— jjPcH c sQ c H eS P c 



2t \ / TljTljj^x 
jx 



t 2 V\l - 2r - 2r 2 ) _ / n m\ (4 ' 3) 

+ £ 2 t / (1 ^ r)2 ES/(son i -s.no), 

where Sj is the spin at site j and the summation i(0) is restricted over the nearest neighbors 
of site 0, only. 

One notices that the first term on the right-hand side when added to P C H C P C yields the 
well-known t-J model, i.e., 

Hu = ~t £ 4 W + ^ £ ( SjSj+x - (4.4) 

jx,a jx 

with c\ a = cJ CT (l — n irJ ). This is expected since the Hubbard Hamiltonian H c reduces to H t j 
in the strong correlation limit. All terms depending on V describe the interaction of the 



magnetic impurity with the strongly correlated electrons. We regroup the different terms in 
and find 

= AE + //kin + H n + H p + H K + H'. (4.5) 
The different contributions are as follows. The term AE describes a constant energy shift: 

where 7=1 + Zt 2 /(U — e/) 2 is again a renormalization factor. The kinetic energy of the 
conduction electrons is given by H kin , i.e., 



//kin — ^ ' Tj [ C,; rr C^j_, rrr ~\- C^j_ vrr C 



J y^j<T^j + Xa ^j + 'jr:' jo 



t f, V 2 (2 + 2r + r 2 ) x 
j ~ 2 { 2e 2 (l + r) 2 ° 3 \ 



(4.7) 



One notices that due to the impurity the electron hopping is reduced between site and 
its nearest neighbors. The interactions between the strongly correlated conduction electrons 
are described by ZZh, which is of Heisenberg type 



5 2 // v 2 \ (4 - 8) 

Again, the interactions between site and its nearest neighbors are reduced due to the 
impurity. H p is a one-particle potential, which describes the attraction of a hole by the 
impurity as well as its repulsion from a nearest neighbor site 

nV 2 

Hp = ' 2| e/ |(l + r) (1 - no) 

7] = 2 + r + -27- - (2 + 7r + 7r 2 + r 3 ) . 

€f(l + r) z v 7 

The attraction of a hole by the impurity is intuitively clear. If a hole is located at site 
a virtual hop of the / electron onto that lattice site creates a singly instead of a doubly 



occupied state. The kinetic energy of the / electron is therefore increased. The last term 
shows that holes are weakly repelled from the neighboring sites i. 

The term ifx is the analogue of the spin-spin interaction in the Kondo Hamiltonian 

2'yV 2 tV 2 (2 + r) 

Hk = JT-SfSo ~ £(s,o + so,), (4.10) 

U-e f Ue f (L+r) . {Q) 

with s j = 1/2 ^Z a p c' ia a a gC Q g. Since site is doubly occupied in the virtual state, the 
dominant prefactor of the first term is 2V 2 /(U — £/) rather than — 2V 2 /ej as in the case of 



un correlated electrons [cf. Eq. (3T)]. The other terms arise from the application of L t and 
indicate a derealization of the interaction. 

Finally, there is a contribution H' to H^ h \ which describes an interaction specific to the 
strong correlation case 

2V 2 t 2 

H' = 77377-3^ E S/ [Si(l - n ) - s (l - n,)] . (4.11) 
ue f {i + r) . {0) 

The first term implies an antiferromagnetic interaction between the impurity and the nearest 
neighbor sites of site provided there is a hole at site 0. Note that in the half-filled case the 
spins of the impurity and those of the nearest neighbor sites i are ferromagnetically aligned. 
The second term can be considered as a correction to the Kondo-type interaction. When an 
electron at site with spin s hops to an empty site i, the / electron can hop virtually onto 
site without creating a double occupancy. Hence, the antiferromagnetic spin exchange is 
decreased. 



B. Half-filled case 



In the case of half filling rij = 1 for all j, and the effective Hamiltonian ( |4.5| ) reduces to 
H c fi = AE + if h + Hk 

*** + t ,-^-*™ (4 , 2) 



2U J 2(17 -e/) 2e f U 2 (l + r 

jx f 
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where Jj is as in Eq. ( f4.8|) . N$ denotes the number of the lattice sites in the substrate. 
The conduction electrons are described by an Heisenberg antiferromagnet with a reduced 
Heisenberg exchange in the neighborhood of the impurity. The / spin couples via an an- 
tiferromagnet spin exchange to the spin at the lattice site which is located nearest to the 
/ orbital. If the (weak) j dependence of Jj is neglected, the Hamiltonian ( |4. reduces to 
that of Ref. |TJ. 

C. Electron doping 

Next, we treat the case of a system with more than half filling (electron doping). We 



proceed in close analogy to Sec. |IV A| , where the case of hole doping was discussed. In 
the limit of large U, empty sites are excluded here. Let P c denote the projector onto the 
configuration space without empty sites and Q c = l — P c . The effective Hamiltonian reduced 
to the space without empty sites reads [cf. Eq. (|4.1|) 1 

= P c H cS Pc + P c H cS Q c — i- — ~ Q c H cS P c , (4.13) 

E (c> — Qc^csQc 

where is the ground-state energy. Only those parts of H e g contribute to the last term 
to order U^ 1 which generate an empty site. Thus, 

# (e) = P c H cS P c - ^jPc(H t + H^)Q c (H t + H®)P C . (4.14) 

Again, we neglect terms of order V 4 , as well as terms involving three lattice sites. Using the 
relation P c rijP c = 1 + n^riji and regrouping the different terms in we find 

#(e) = ^(e) + jj® + R (c) + H {c) + R {e) ? ^ 15) 

where 
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A£ (e) = e/ - 2( J^ e/) + W - N a ) 

rr(e) _ j-(e) M , , p t p 

kin / y ■ i j V ^ja^j+xa ~ ^j+xa^ja 



JX,<7 



j°j+x 
jx 

-Et( 1_ n 3l n n)( l ~ n j+x ^n j+xl ) 



(4.16) 



4 

JX 



2(U-e f ) 

r(e) 



H*> - 271/2 Ss 1 tV * S Vfs +s 1 



N$ is the number of lattice sites in the substrate, N c is the number of conduction electrons, 
so that there are N c — N s doubly occupied sites, which contribute with an energy U(N C — N s ) 
to the energy shift AE^ C \ 

As in the case of hole doping, the second and third term in Eq. ( |4.15|) , + define 
a t- J model describing the conduction electrons with parameters, which are modified in the 
neighborhood of the impurity: 



T (e) - - M V 5 

lj ~ "2 \ 1 " 2{U-e f )^\ 



J i~lT\ l -\e f \{U-e f ) 5 * 



(4.17) 



In the case of electron doping, double occupancies hop rather than holes: cJ CT = c^n^. Since 
the / electron can gain kinetic energy only if site is singly occupied, double occupancies are 
repelled from site 0. This is described by the "potential", H^°\ The Kondo-type interaction, 
, between the / spin and the spin of the conduction electrons consists of two parts. The 
local part is the same as in the case of hole doping, cf. Eq. ( |4.10| ). The prefactor of the 
second term is modified and we introduced s 0i = 1/2 J2 a j3cl a cr a/3 c 0f3 . However, in the case 
of electron doping there are no interactions of the / spin with the conduction electrons 
corresponding to H', cf. Eq. ( 4.11 ). 



Finally, we mention that in the half-filled case n^riji = 0, and the Hamiltonian ( (4.15 ) 
reduces to that of Eq. ( |4.12|) . 
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V. SUMMARY 



In this paper, we considered a magnetic impurity which interacts weakly by hybridization 
with a system of strongly correlated electrons. This contrasts the situation treated in the 
Anderson impurity model where the impurity interacts with free conduction electrons. We 
described the strong correlations among the conduction electrons by a Hubbard Hamiltonian 
and performed a canonical transformation to eliminate the charge degrees of freedom at the 
impurity site. For vanishing Hubbard repulsion U, this procedure was shown to be equivalent 
to the conventional Schrieffer- Wolff transformation and yields the Kondo Hamiltonian. 

Of particular interest is the limit U >> t. Here, we additionally reduced the charge degrees 
of freedom of the conduction electrons by a degenerate perturbation expansion which, in the 
absence of the magnetic impurity, leads to the t-J model. The corrections we find are due 
to the impurity. In the case of half filling, the conduction electrons are described by a 
Heisenberg Hamiltonian with an exchange coupling constant, which is slightly modified in 
the neighborhood of the impurity. Furthermore, there is an exchange interaction between 
the impurity spin and the spin of the conduction electron at the lattice site closest to the 
impurity, i.e., site 0. The coupling constant depends on U and is, therefore, different than 
in the usual Kondo Hamiltonian. 

When the system is doped new terms arise. Since we strictly forbid double occupancies 
of the impurity site the electron-hole symmetry is broken and we have to treat the cases of 
hole and electron doping separately. If the system is doped with electrons, doubly occupied 
sites are repelled by the impurity since otherwise the exchange becomes impossible. 

When we dope the system with holes, the conduction electrons move like in a t-J model. 
The holes are attracted towards the impurity, which reflects the gain in kinetic energy if the 
/ electron hops without creating a doubly occupied site. Additionally, we find a new type 
of interaction: Provided a hole is located at the lattice site there is an antiferromagnetic 
exchange of the impurity spin with the conduction electron spin at the nearest neighbor 
sites of 0. In the undoped system these spins are ferromagnetically aligned. 
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APPENDIX A: EFFECTIVE HAMILTONIAN 



In the following, we derive the effective Hamiltonian ( 2.12| ) which is obtained from the 
starting Hamiltonian (|2.4| ) after performing the canonical transformation (|2.5|) that elim- 
inates the charge degrees of freedom of the / orbital and projecting onto the space with 
rif — 1. We begin by listing the terms to lowest order in the expansion ( j2.1U| ) of l/(ej — L c ). 
The zeroth-order term [y = 0) reads 

1 + 1 / 1 



For v = 1 one obtains 



e f -L v e f -L 



u 



'Oct 



-Lr 



-y 

-/ i(0) 



'Oct 



ill 



1 + r 



(Al) 



1 - 



1 + r 



1 + r 
1 



1 + r 



n 05 a 



1 + r 



n i<T I C \a C 0a C 0a + 



l + r) 1 



(A2) 



and for z/ = 2 



-f i(0) 



' (1 + 2r) 2 

1 " 7TW nos + 



'1 _|_ r )2 niCT y Co<J (1 + r )2 V WsU ° sWct ^ctWctWct, 

The summations z(0) are restricted over nearest neighbor sites i of site and in Eq. ( |A3| ) we 
have neglected contributions from next-nearest neighboring sites. When the Eqs. (|A"1|- [A3| ) 
are set into Eq. (|2.9|) one obtains H c r in the form 



(A3) 



H eS = H c + e f + H®+HW+HW. 



(A4) 



The different contributions are 



V 2 



H^ = V -- V ^ + r) n Q + 

6/ 26/(1 + r) e^l + r) 



^ot^oi + 



2\/ 2 



t/(l + r) 



S/s , 



(A5) 
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and 



H (2) 



tv 2 



4e 2 ,(l+r) 



Z ( c l c Oa + c L c ia) (2 + 2r + r 2 - (2 + r)n ^ - (3 + r)n is + 3^?^ 



-'icr L '0(T 1 ^Oa^ia I l i i ' 1 ' J '"V(7 V." 1 ' ) '"la i <-""U<t"'ict , 

+ o an Z Z ( c L c o<7 + cLc iCT ) (r(r + 2) + (1 - r)n i(f - rn 0ff - n 05 n i(f ) 

/e^i + rj i(0)o . 



t\/ 2 



26}(l + r) 2 i(0) 



Z (^/ s io + S f 4i) ( r ( r + 2) + (1 - r)n iT - rn | - n oi n iT ) 



2e} ( 



tv 

TTXZW Z ( 5 / % + 5 / s io) ( r ( r + 2) + (1 - r> a - rn 0T - n 0T n a ) 
l i + r J i(0) 



(A6) 



ZV 2 t 



2+2 



VH 



2+2 



2Ue){l + r) 



c / "~ "-/V- 1 ' / i(0) 

Z\/ 2 t 2 (2 + 7r + 7r 2 + r 3 ^ 



Z ( c 1t c 1| c o| c ot + c OT c o| c i| c iT 



VH 



2+2 



2e 3 (l + 



-n 



Zl/ 2 t 2 (l + 2rf 



\/ 2 t : 



2+2 



z /V x -rrj ^c /v ^ , , , . (0) 



2[/e 2 f (l + r) 2 



Z n * H ^ 3 ' n 0T^0| 



e 3 f (l + r) 



£/e 2 (l + r) 2 ^ 



(0) 



Z ( s o s * - 



4 



V 2 t 2 _ Q /2(1 + r-r 2 ) \ ^ 2 t 2 _ a . 

Z S / s o V , , ^ -n l ) + — m ■ , 2 Z S / S » (2 -no) 



i(0) 



1 + r 



i(0) 



(A7) 



The Hamiltonian if e g is discussed in several limiting cases in Sees. |T| and [TV 
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